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Abstract. The three dimensional Gribov-Zwanziger Lagrangian is analysed at one and two 
loops. Specifically, the two loop gap equation is evaluated and the Gribov mass is expressed in 
terms of the coupling constant. The one loop corrections to the propagators of all the fields are 
determined. It is shown that when the gap equation is satisfied the Faddeev-Popov ghost and 
both Bose and Grassmann localizing ghosts all enhance in the infrared limit at one loop. This 
verifies that the Kugo-Ojima confinement criterion holds to this order and we also show that 
both Grassmann ghosts are enhanced at two loops. For the Bose ghost we determine the full 
form of the propagator in the zero momentum limit for both the transverse and longitudinal 
pieces and confirm Zwanziger's recent general analysis for the low energy behaviour. We pro- 
vide an alternative but equivalent version of the horizon condition expressing it as the vacuum 
expectation value of an operator involving only the localizing Bose ghost field. The one loop 
static potential is also determined. 
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1 Introduction. 



The Gribov-Zwanziger Lagrangian is a formulation of the Landau gauge fixed Yang-Mills theories 
where the Gribov problem is incorporated in a localized way, [JJ [2j O 0] . This problem, [5], 
essentially relates to the difficulties in fixing a gauge globally for gauge theories with a non- 
abelian symmetry. In his seminal work, [5], Gribov demonstrated that globally different gauge 
configurations could satisfy the same gauge condition thereby introducing an ambiguity into the 
gauge fixing procedure. Such Gribov copies do not affect the local gauge fixing in Yang-Mills 
theories and hence the ultraviolet structure of such theories does not encounter gauge fixing 
difficulties. By contrast, the problem relates to global issues and hence the infrared regime 
of the theory. For non-abelian gauge theories, Gribov indicated that such gauge copies could 
be entangled with the problem of confinement, [5], which is sometimes referred to as infrared 
slavery. One consequence of the analysis of [5j is to overcome the copy problem in the main 
by restricting the path integral to a specific region of configuration space. This region, known 
as the Gribov region, denoted by SI and containing the origin, is defined by the locus of points 
where the Faddeev-Popov operator is positive, [5]. Geometrical aspects of the region and their 
consequences have been explored in [6j [8j [91 \TU[ [TT] . As an aside we note that such a restriction 
does not lead to unambiguous gauge configurations. Instead the Gribov region has a subregion 
called the Fundamental Modular Region, denoted by A, where the gauge is fixed uniquely 
globally. However, it has been argued in [TT] that Green's functions defined over A and are 
equivalent. To incorporate the path integral restriction to Q Gribov modified the Yang-Mills 
action to include a non-local operator which in effect cut off the domain of integration, [5]. 
The presence of such an operator, referred to as the horizon or no pole condition, introduces 
an arbitrary mass scale, 7, which is known as the Gribov mass. However, it is not a new 
parameter of the theory but satisfies a gap equation defined by the defining horizon condition 
and is a function of the coupling constant, [5]. The presence of this non-local operator and the 
Gribov mass alters the structure of the propagators of the theory. For instance, the gluon has 
a propagator which vanishes at zero momentum and depends on the Gribov mass. Though it 
has a non-fundamental form with the gluon being in effect massless but with a non-zero width, 
[5]. This may appear to be contradictory but the gluon is not a fundamental field in itself as it 
is confined. A second feature is that as a consequence of the gap equation for 7, the Faddeev- 
Popov ghost propagator has an enhanced or dipole behaviour in the zero momentum limit. The 
latter property was later encapsulated in the Kugo-Ojima confinement criterion, |12|. I13J . for the 
Landau gauge. Indeed this condition has been re-examined in the Gribov-Zwanziger context in 

The relevance of the Gribov-Zwanziger Lagrangian to the Gribov path integral restriction 
rests primarily in the reformulation of Gribov's Lagrangian in a localized way, [JJ [21 [3j Hj • Gri- 
bov's analysis was in a semi-classical approximation but one cannot perform high order loop 
computations using a Lagrangian with a non-local operator. Instead Zwanziger managed to lo- 
calize the non-locality to produce a local renormalizable Lagrangian for the Landau gauge. The 
renormalizability has been established in several articles, I15|. [16] . The localization introduces 
several extra fields which are known as ghosts. However, one set, {(/>^ fe , <^f}, have Bose statistics 
whilst their partners, {00^ , <D° 6 } , are fermionic. The latter are crucial in maintaining the estab- 
lished ultraviolet properties of the theory such as asymptotic freedom and ensure the one and 
higher loop MS j3- function, [TTJ [TBI 032 [2Q1 EEj , is unaltered. This is important since the Gribov 
operator in some sense relates to the infrared structure of the theory and its presence therefore 
ought not to upset the ultraviolet structure where indeed there is no gauge fixing ambiguity. 
One consequence of the localization is that one can carry out explicit loop computations. In 
[21 U] the one loop gap equation satisfied by 7 in the original Gribov action was reproduced. 
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This has been extended to two loops in [22J in the MS scheme and led to the check that the 
Kugo-Ojima confinement criterion holds to two loops explicitly. More recently the one loop 
static potential for heavy quarks was computed in the Gribov-Zwanziger context in [23] as well 
as the full one loop corrections to all the propagators of the fields. It transpires that the latter 
have been important for verifying a recent non-perturbative analysis of the Gribov-Zwanziger 
framework by Zwanziger, |24| . 

Whilst the appearance of dipole propagators for the Faddeev-Popov and a;" ghost propaga- 
tors is in keeping with the Kugo-Ojima criterion, [121 [T3l 114] . such fields cannot play a role in 
the actual confinement of heavy quarks, for instance. This is purely due to their statistics and 
the lack of a (direct) coupling to quarks. Instead one would require an enhanced field with Bose 
statistics. Clearly the gluon cannot be that field due to its infrared suppression. However, using 
symmetry arguments which are valid to all orders, Zwanziger has argued that certain colour 
components of the imaginary part of the Bose ghost field, </>^ fe , are enhanced, [21]. Indeed this 
structure has been confirmed at one loop in [25J. There the explicit enhancement was shown 
for the transverse component. The longitudinal piece was not considered since it clearly could 
not play a role in the exchange particle for the static potential considered there. However, it 
has been analysed subsequently (and referred to in passing in [21]) and will be reported on in 
full in this article as part of a larger calculation. More specifically given the elegance of the 
Gribov-Zwanziger construction and its potential for being a working Lagrangian incorporating 
confinement, it is the main aim of this article to record in one place the full analysis for the 
three dimensional theory. In the series of papers, |22} [23"1 125j . the two loop MS gap equation, 
static potential and all the propagator corrections were all given for four dimensions. However, 
if one is to have a full understanding of that case, it should also be the case that there are 
parallels in the lower dimensional version. Indeed the three dimensional theory has several in- 
teresting features deserving study in their own right. For instance, the Gribov mass plays the 
role of a natural infrared regulator in this superrenormalizable quantum field theory. Moreover, 
the ultraviolet finiteness means that the gap equation simply relates the Gribov mass to the 
(dimensionful) coupling constant. Therefore, we will provide all the quantities for the three 
dimensional Gribov-Zwanziger Lagrangian that have been computed in four dimensions to the 
same loop order. 

It would be remiss not to discuss the relation of the current Gribov-Zwanziger scenario with 
that found on the lattice for quantities such as the gluon and Faddeev-Popov ghost propagators. 
The present point of view is that the propagators are not respectively suppressed or enhanced at 
zero momentum. Instead the gluon propagator freezes to a non-zero value whilst there is clearly 
no dipole behaviour for the ghost. The evidence for this has been provided over a number of 
years by several lattice collaborations, [261 123 [23 [291 USB Ell E2] • This decoupling scenario, [33], 
is in contrast to the conformal or scaling situation of gluon propagator suppression and ghost 
enhancement of the original Gribov set-up, [5]. Whilst it is possible to model the decoupling 
situation by a condensate argument based on the Gribov-Zwanziger Lagrangian, [Ml [35], it is 
clear that the Kugo-Ojima confinement criterion cannot be satisfied. Moreover, there cannot 
be any enhanced propagators with either Bose or fermion statistics. However, one can argue 
that the positivity violating gluon propagator which the decoupling solution has, is sufficient to 
ensure a confining theory. Though a condensate explanation based on a perturbative vacuum 
would need to be extended to incorporate non-perturbative aspects of the vacuum. Irrespective 
of this we believe that the debate has not been fully resolved and that to understand the theory 
one ought at the very least to have as much analysis available as is calculationally possible and 
in spacetime dimensions other than just four. 

The article is organised as follows. We review the main features of the Gribov-Zwanziger 
Lagrangian in section two which are required for our three dimensional analysis. The two loop 
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correction to the gap equation satisfied by 7 is given in the next section with an estimate for 
the non-zero one loop value of the renormalization group invariant strong coupling constant 
at zero momentum. Section four is devoted to the calculation of the one loop static potential 
of heavy quarks. Given that there is currently interest in the behaviour of the Bose ghost at 
zero momentum we give the formal one loop propagator corrections in section five. Whilst the 
transverse part was considered in |25| . we concentrate on the longitudinal piece and extract 
the Landau gauge behaviour in accord with Zwanziger's analysis of [24]. The explicit one loop 
structure of the gluon and propagator sector for both three and four dimensions and for 
arbitrary colour group are recorded in section six. We give our conclusions in section seven. 
There are three appendices. The first two record the explicit one loop corrections to all the 
2-point functions for the transverse and longitudinal sectors respectively. The final appendix 
provides the complete structure of the one loop form factors appearing in the propagator of the 
real part of (f> ab given the most general possible SU(N C ) colour structure of the corresponding 
2-point function. 



2 Formalism. 

In this section we recall the basic formalism for the Gribov-Zwanziger Lagrangian, [U [21 [3j H]. 
First, the canonical QCD Lagrangian with a linear covariant gauge fixing term is given by 

L QCD = _ \ Q a na^v _ J_( d H A a^2 _ giQUn c * + ^ (2.1) 

4 M 2a p 

where Cr" is the field strength for the gauge potential A^, the covariant derivative, D^, is 
defined by 

D^c a = cV a - gf ahc A\c c 
D^ u = d^ u + igTfjA^ (2.2) 

g is the coupling constant, f abc are the colour group structure constants with generators T a . The 
various indices are restricted to the ranges 1 < a < Na, 1 < I < Np and 1 < i < Nf where Nf 
and Na are the respective dimensions of the fundamental and adjoint representations and Nf is 
the number of massless quarks. Whilst the focus is primarily on the Yang-Mills Lagrangian we 
have incorporated massless quarks partly for completeness and also as an internal aid in checking 
some computations. The gauge parameter a is included in order to assist with determining the 
propagators of the theory. However, we will work in the Landau gauge throughout, where a = 0, 
which is assumed unless it is required at intermediate stages of computing one loop propagator 
corrections as will be the case in a later section. The Lagrangian (12. ip is the usual starting 
point for high energy computations and one does not need to be concerned with the fact that 
the gauge is not fixed uniquely globally. In the ultraviolet regime Gribov copies do not alter 
physical predictions. However, to handle the ambiguity problem the path integral restriction 
equates to modifying the action by the no pole condition or equivalently the horizon condition. 
The boundaries of the Gribov regions are given by the zeroes of the Faddeev-Popov operator 
d^D®. So the interior of the first Gribov region is that set of points where the inverted Faddeev- 
Popov operator is finite. Whilst the original arguments of [5] were based on a semi-classical 
approach this now equates to extending (|2.ip to the Lagrangian 

L Grib = L QCD O^t 1 _ dNA^ 

2 d v D v 2g y 1 



where 7 is the Gribov mass parameter. Originally the non-local operator of (|2.3p was only 
approximated by the Laplacian, [5], but this was later extended and made more concrete by 
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Zwanziger in [lj. The parameter 7 is not an independent quantity in the Gribov theory. Instead 
it is a function of the coupling constant and the relation between the two is defined by the 
horizon condition. For (I2.3p this is, [HE], 

W^A-M) _ (2.4) 

where Ca is given by 

jacdfbcd = CA§ ab ( 25 ) 

and d is the spacetime dimension. If one could handle the non-locality when calculating the 
vacuum expectation value of (|2.4p then the gap equation satisfied by 7 would emerge. In four 
dimensions 7 is expressed as a non-analytic function of the coupling constant. It is important to 
stress that one cannot treat 7 as an independent parameter of the theory. The non-local theory 
cannot be regarded as a gauge theory, in the Landau gauge, unless 7 satisfies the gap equation, 

The key to resolving the calculational obstacle represented by the non-local term of (|2 . 3[) was 
provided by Zwanziger in a series of interrelated articles, [HE1E1II1E1IZ1E1E]. By considering the 
properties of the Gribov region in the Landau gauge the non-local Lagrangian was transformed 
into a local Lagrangian which involved extra spin-1 fields. These localizing ghosts, 4> ab , 4> ab ', uj® b 
and O^, where the first pair are bosonic and the latter Grassmannian, are additional to the 
gauge potential and the Faddeev-Popov ghosts. Their presence does not alter the ultraviolet 
properties of the theory since, for instance, asymptotic freedom still holds in four dimensions. 
Instead they become effective in the infrared limit as one approaches the Gribov boundary. More 
specifically, the localized version of the Gribov Lagrangian is [H El El HI El El EE] 

L GZ = L QCD + ± f f*Hff'{p vPll ) ab + l -p ab ^(D u ^) ab ~ l -e b ^d u {D vPtl ) ab 

- -l=gr b ^^(D u c) b e cfl - il 2 f abc A a ^ c - (2.6) 

where there is a mixed 2-point term involving the gluon. We have chosen to follow the current 
convention and use the real and imaginary parts of the Bose ghosts rather than the complex 
versions, [24, 36J. This is because in four dimensions the behaviour of the propagators of each 



formulation. We take as the real and imaginary parts 



component is significantly different and is difficult to extract cleanly in the original (j) ab and r; , 



"'' " 7^ 0? + O , 4f = ^= (p? - O ■ (2.7) 



A* 



(For comparison p ab and are respectively the U® b and V^ b fields of [MJ EE].) Although we 
have omitted the a dependent term since (|2.6p corresponds to the Landau gauge, one requires 
that term to safely derive all the propagators. As we will be considering the infrared properties 
of the one loop corrections to the transverse and longitudinal parts of the propagators, we record 
for completeness the form of the intermediate propagators prior to taking the a — > limit. These 
are 



5 ab n 2 a5 ab p 2 



be if abc j 2 iaf abc ^ 2 

(A«(p)C(-p)) = [(p2) 2 + C - A7 4] i V(P) + [ (j)2) 2 + aQ7 4 ] V(ri 



5 



{A a Jp)p b v c { 
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-P)> 




-P)) 







= 

= «(p)a)^(-p)) = - -jp-Vp, (2.8) 

P^(p) = - ^ , V(p) = ^ (2.9) 

are the respective transverse and longitudinal projectors. Consequently, in our Landau gauge 
calculations we will use 



where 



(Al(p)Al(-p)) = - —^JL^P^fr) 



5 ab p 2 



(A«(p)p b u c (-p)) = 

xacxbd sabe fcde^.4 

(c b (P)c d (-P)) = - + p2[( ; 2) / + j A74] ^(P) 

(Cf(p)pt d (-P)) = o 

xacxbd 

(P»(p)pl d (-P)) = «(pR d (-p)> = ~ —r,, u (2.10) 



as our propagators. The derivation of (|2.8|) and (|2.1Up is complicated by the mixed term of 
but was discussed at length in [25] . Though we note that for the real and imaginary Bose ghost 
there is a clean split of the propagators with the real part propagator having a similar form to 
that of the associated fermionic localizing ghost. When we examine the propagator corrections in 
the infrared this feature will be preserved in keeping with Zwanziger's recent general arguments, 



3 Gap equation. 

In this section we derive the two loop gap equation satisfied by the Gribov mass in three 
dimensions. This computation is similar, for example, to that in four dimensions in terms of the 
Feynman diagrams to be computed. The method is to evaluate the horizon condition (|2.4ft but 
not in the non-local version of the theory. Instead we consider the equivalent definition of the 
condition in the localized theory, (12. 6h . In terms of the real Bose ghost fields this is 

r hc (A^{ X )il\x)) = (3.1) 

where the relation between the A ^ and ^ is established via the equation of motion 

A% = - -^^{d^D^t . (3.2) 

Hence (|3.ip clearly equates to f|2.4[) using (|3.2p . Using this version of the gap equation one 
evaluates the Feynman diagrams of the vacuum expectation value to two loops. At leading order 
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there is one Feynman graph and at two loops there are nineteen diagrams to evaluate. These 
are generated using the Qgraf package, [37], and then converted into Form input notation 
where Form, |38j, is the symbolic manipulation language used to handle the associated algebra 
with the computation. We follow the standard procedure of breaking the one and two loop 
Feynman diagrams up into a sum of master vacuum bubble integrals by using tensor reduction 
and then substituting their explicit values. For three dimensions, all master massive vacuum 
bubble topologies to three loops have been computed in [39] for all possible independent masses. 
It is relatively straightforward to extract the integrals required and include them in the Form 
routines. However, one needs to be careful in the Gribov-Zwanziger case where the Gribov 
propagator is not the standard massive propagator. One first has to apply partial fractions 
using, for example, 



1 



[(p 2 ) 2 + Cat 4 ] 2zV^7 2 



1 



\p 2 



(3.3) 



to obtain factors within the Feynman integrals of more standard form. However, each involves 
an imaginary mass corresponding to massless unstable fields. The main issue, though, is in 
utilizing the master integrals of [39] which we illustrate with the simple one loop integral. From 
[39] 



d d k 



1 



m 
4tt 



1 + 



2 + 21n 



1 

2m 



+ 0(6 2 



(3.4) 



(2n) d [k 2 + m 2 } 

where in dimensional regularization d = 3 — 2e. The argument of the logarithm will be rendered 
dimensionless when the mass scale, /i, which is required to retain a dimensionless coupling 
constant in (f-dimensions is included in the overall computation. In four dimensions the overall 
factor of the evaluation would be dimension two but in three dimensions the dimensionality 
reduces by one unit. Thus for the Gribov case one would require the square root of the width. 
For each of ±iv / C ! A7 2 there are two possibilities resulting in four different underlying masses. 
However, to exclude any potential ambiguity in the overall final expression for the gap equation, 
which must be real and not complex, within our Form routines we have formally extended (|3.4p 
to 



d d k 



1 



/ C/i7 2 



(2vr) d [k 2 + iy/C^-f 2 } 



4tt 



1 + 



2 + 2 In 



1 



e + 0(e 2 



(3.5) 



and its complex conjugate. Then in constructing our final overall gap equation we merely use 
the unambiguous and elementary identifications 



(3.6) 



If these objects appear instead in a ratio then one first rationalizes before using these trivial 
identities. For completeness we note that the basic two loop sunset topology in three dimensions 
with three distinct masses mi, and is, [39j, 



1 



ki [k 2 +m 2 ][l 2 +m 2 ][{k-l) 2 + 



16vr 2 



1 1 . 

1 h In 

4e 2 



1 



[mi + m 2 + m 3 ] 



+ 0(e) 



(3.7) 



which is the only non-trivial topology at two loops. The other main topology is the product of 
two one loop vacuum bubbles. 



The one loop integral, (I3.4p . illustrates one feature of the three dimensional Gribov-Zwanziger 
Lagrangian which differs from the four dimensional case and that is that to two loops the theory 
is ultraviolet finite. Though master integrals, such as (|3.7p . can be divergent. Ordinarily in 
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three dimensional Yang-Mills theory with a massless gluon one has to be aware that the theory 
is potentially infrared pathological. However, in (|2.6p the presence of the Gribov mass whilst 
not corresponding to a non-zero mass does act as an infrared regulator. This is relevant at 
two loops for the computation of the gap equation since at one loop only (|3.4p is required. At 
two loops within our routines we have been careful in noting the potentially infrared divergent 
master graphs and checked that they actually cancel among themselves in the overall sum of all 
the contributing integrals for a diagram. These arise, for instance, when there are propagators 
of the form l/(/c 2 ) 2 which occur in the propagator due to the transverse projector. Whilst 
it is clear that overall such a factor has a tensor structure rendering the potential double pole 
as unproblematic, within the algebraic rearrangements to produce the master integrals one may 
have similar problematic integrals which are only cancelled, for example, by using integration 
by parts. Although the finiteness of (12. 6h may appear beneficial from a computational point of 
view it is actually a disadvantage. In a renormalizable ultraviolet divergent theory, such as the 
four dimensional version of (I2.6p . the renormalization constants satisfy Slavnov- Taylor identities, 
[H fl"5"| [16] . These then serve as important checks on setting up the Feynman diagrams and the 
actual explicit master integral evaluations. In the three dimensional case we do not have these 
additional important internal checks. However, the main parts of the Form code are the same 
as the four dimensional work and we merely use these again as they have been checked. This 
leaves us with minimizing the potential source of errors to be that of substituting for the master 
integrals correctly. Of course, for the gap equation we must have a real expression ultimately 
since 7 is a real parameter. It is not an independent quantity since it is defined by the horizon 
condition and will be a function of the coupling constant which is real. In three dimensions, of 
course, the coupling constant is dimensionful and hence the gap equation will relate quantities 
to ensure that overall there is only one independent dimensionful parameter in (|2,6p . 



Given these considerations the two loop gap equation for 7 in (|2.6p is then 



3 = VlcTg 2 + 

4 16-7T7 
.6 



917vr 17 545 

H - H — tan 



262144 98304 131072 



C A — T F Nf 

256 1 



C 1/2 n 4 

7T 2 7 2 

+ 0(5°) (3-8) 

for Nf massless quarks where Ca is the usual adjoint Casimir. The arctangent derives from 
the four complex masses defining the Gribov propagators and the form of the finite part of 



(|3.7p when, for instance, there are two propagators giving a mass J iyj CaI 2 and one with 



iy / C J 47 2 . This produces a term ln(4 + 3i) and its conjugate for the conjugate integral and 
within the overall computation it is the imaginary part which is translated into the final gap 
equation. As a note on our conventions each appearance of 7 is always with one factor of C 1 / 4 



so that the peculiar appearance of these factors and powers in (|3.8p is actually consistent with 
the presence of Ca which is what ordinarily appears from the group theory in the one loop term. 
As the three dimensional theory is ultraviolet finite then both the coupling constant and 7 do 
not run. Moreover, there are no logarithms involving 7 and the renormalization scale \x as there 
is in the four dimensional gap equation. In [22] the four dimensional gap equation was solved 
in order to write 7 as an explicit function of the coupling constant producing an explicitly non- 
perturbative function. For (j3.8|) we can also relate these parameters. The way we have chosen 
to do this is to simply solve (|3.8p as a quadratic equation. As noted in [25] there does not 
appear to be a unique way of solving the gap equation. Choosing to solve as a quadratic here 
is straightforward but if the explicit three loop or higher gap equation was known then it is 
not clear whether a numerical solution could be extracted in those cases. Despite these caveats 
we have solved (|3.8p numerically for both SU(2) and SU(3) for a variety of values of Nf and 
recorded the results in Tables 1 and 2. More specifically we have introduced the dimensionless 



S 



variable A n defined by 



A, 



(3.9) 



where n reflects the loop order. Each table provides the relation of 7 to the coupling constant 
and vice versa. Clearly from the tables using this method of solution the convergence does not 
appear to be very good. Although SU(3) is better than for SU(2). It would be interesting to see 
if the three loop corrections improved the convergence but one glance at the explicit expression 
for the master three loop integral for the Benz vacuum bubble topology of [39] would indicate 
how tediously complicated such a calculation would be. 



N f 


Ai 


A 2 


Ar 1 


A^ 





1.60660 


0.60389 


0.62243 


1.65593 


2 


1.60660 


0.70962 


0.62243 


1.40920 


3 


1.60660 


0.79534 


0.62243 


1.25732 


4 


1.60660 


0.93636 


0.62243 


1.06797 



Table 1. Numerical values for relation between g and 7 for SU(2). 





Ai 


A 2 


a: 1 


A 2 





0.70711 


0.40260 


1.41421 


2.48385 


2 


0.70711 


0.44502 


1.41421 


2.24709 


3 


0.70711 


0.47308 


1.41421 


2.11381 


4 


0.70711 


0.50851 


1.41421 


1.96653 



Table 2. Numerical values for relation between g and 7 for SU(3). 



Given we have obtained a relation between the mass parameter, 7, with the coupling constant 
from a two loop gap equation in three dimensions, it is worth noting related work. One interest in 
three dimensional Yang-Mills theory resides in the fact that it is relevant to the four dimensional 
finite temperature theory. In this situation it is believed that a non-zero magnetic mass is 
generated dynamically non-perturbatively. Such a magnetic mass can be accessed from gap 
equations. For instance, there are a variety of one loop results available, [3UJ HU |32j |33] , where 
[41 j. [42] used a non-local mass operator but which was not of the Gribov form. These ideas 
were extended to two loops in [33]. There an estimate of the magnetic mass, m m , was quoted 
as m m 0.34g 2 for SU{2). In our case for the case with no quarks C]/ 4 7 = 0.132g 2 where the 
group factor is unevaluated and included with 7 because of our conventions. The discrepancy in 
values here ought not to be taken seriously though as in the former case the method of attack 
is to have a canonical massive gluon propagator. By contrast we are considering a Gribov style 
of propagator where the mass is zero but the width is not. 

Next we use the one loop gap equation to determine the leading order value of an effective 
coupling constant which has been shown to freeze to a finite value, [24]. From the gauge 
potential and Faddeev-Popov ghost propagator form factors one can define a renormalization 
group invariant object which behaves as the coupling constant at high energy. This is primarily 
due to the fact that the gluon ghost vertex does not undergo any renormalization due to a 
Slavnov- Taylor identity, [45]. The source of the zero momentum value being non-zero derives 
from the momentum dependence of the form factors when the gap equation for 7 is set. More 



9 



specifically, if we define the gluon propagator as 



(A«(p)Al(-p)) = ~ 5 ab ^^P, u (p) (3.10) 



where Da(p 2 ) is the form factor and 

(c a (p)c b (-p)) = 5 o6 ^P (3.11) 

for the Faddeev- Popov ghost then the renormalization group invariant effective coupling constant 
is defined by 

af{p 2 ) = a s ^)D A (p 2 )(D c (p 2 )) 2 . (3.12) 

We have computed the one loop corrections to both Da(p 2 ) and D c (p 2 ) for (|2.6p in three di- 
mensions and recorded the explicit functions for each in Appendices A and B. That for the 
Faddeev-Popov ghost is equivalent to due to the similarities between the real Bose ghost 
and fields. This follows from the consequences of the underlying Slavnov- Taylor identities 
for (|2.6p which have been discussed primarily in the context of the four dimensional theory but 
which also are valid in the three dimensional case, [H 115} ITS]. Therefore since the gluon propa- 
gator vanishes as 0{p 2 ) as p 2 — > 0, meaning Da{p 2 ) is 0((p 2 ) 2 ), and D c {p 2 ) is also 0(p 2 ) in the 
same limit one is left with a finite answer for the effective coupling constant at zero momentum. 
Specifically we have 

of (0) = ^cTl ■ (3-13) 

We recall that since we are in three dimensions the coupling constant carries a dimension which 
is why 7 appears on the right hand side. However, 7 is not independent and satisfies the gap 
equation being reexpressed as a function of the dimensionful coupling constant. Although this 
is a leading one loop calculation and therefore qualitative, the three dimensional set-up may 
be more useful in exploring the Gribov-Zwanziger scenario further. For instance, if one could 
obtain a non-zero estimate for a frozen effective coupling constant, then this would essentially 
fix the parameters of the theory provided the gap equation was known sufficiently accurately. 
The latter would be essential if, for instance, one wanted to extract a reliable magnetic mass 
estimate. 

We close this section with an indication of an alternative way of computing the gap equation. 
In (|2.3p the definition of the horizon condition, (|2.4p . involves the non-local operator which 
cannot be determined without localization. Reformulating (|2.3p in terms of localized fields 
produces a local version of the horizon definition, (|3.ip . by virtue of ()3.2|) . Given this we can 
reformulate (|3.ip again by eliminating A® within the vacuum expectation value to produce an 
expectation involving only fields at leading order. In other words 

f abp f cdp (e b »{x){d»D v ^) cd {x)) = - d ° A ^ (3.14) 



should also be equivalent to the horizon definition and produce the same gap equation as ([3. 81) . 
We emphasise that this gap equation is not the vacuum expectation value of the kinetic 
term due to the presence of the structure constants. So there is no parallel definition for p ab . 
Given this reasoning we have evaluated the one one loop and eighteen two loop vacuum bubble 
graphs contributing to (|3. 14[) . This uses the same basic one and two loop master integrals as 
that for (|3.ip and it is satisfying to record that (|3. 14[) does indeed reproduce f|3.8j) . Moreover, 
we have also checked that the two loop MS gap equation of [22] in the four dimensional theory 
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is also recovered with the definition (I3.14p . This is more involved than the three dimensional 
case as one has to correctly take account of the ultraviolet divergences. So we can summarize 
the different definitions of the gap equation in the unifying equivalences 

(^ (x) s4/ aM(x) ) = -c^r bc {A a "{x)e, c {x)) 

= -^^(e b ^x)(d"D^) cd (x)) = i^L. (3.15) 

In the last definition, which involves two terms when the covariant derivative is written explicitly, 
we have used (|3.2p to redefine the gluon field within the vacuum expectation value. However, 
there is in principle no reason why one cannot repeat the substitution of (|3.2p in the covariant 
derivative of (|3.14p . This would produce a vacuum expectation value involving three terms and 
equate to a perturbative expansion where the final term will involve a gluon via the appearance 
of a new covariant derivative. Iterating this procedure one can replace the final horizon definition 
by an infinite series of terms involving only fields in a perturbative expansion. For instance, 
the first few terms would be 

dCANAlA = ft d (d v e b »[d„tf - -^ft s (d°d a e„ s )tl d 



9 2 



"' 2 -Jt frs fl qmn d° Wd P C n ) &] tl d 



ch 4 

+ 0( 5 3 ) ) (3.16) 



where f^ bcd = j ab Pj cd P t Q simplify notation and we have integrated by parts on the ordinary 
derivative. With this version of the horizon definition we have evaluated the one one loop and 
eighteen two loop graphs contributing to (|3.16p in three and four dimensions and found that 
the respective two loop expressions for the gap equation are reproduced exactly. Within the 
context of the vacuum expectation value definition of the horizon condition one might regard 
the perturbative expansion of (13. 16f> as an infinite series representation of the original non-local 
definition of Gribov, (|2.4p . In other words in (|3.15p the first and last vacuum expectation values 
are a field theoretic type of geometric series with regarded as a pseudo-dual field to the 
gluon. 

Whilst these equivalences between the different formulations of the horizon equation are 
novel and suggest a type of duality between the A^ and fields due to (I3.2p . one must be 
careful when it should be set. For instance, it is tempting to replace all appearances of in 
the localized Lagrangian, (12.6p . in the hope of producing some sort of effective low energy field 



theory where the dominant field is £ t ® b . However, if one naively does this then the propagators 
of the theory have no relation to that of ^ in (I2.10p . For instance, eliminating the gluon 
completely in this naive approach would produce the non-standard propagator, for non-zero a, 

sacxbd fabe fcde 

<£»€ d (-p)> eff = - -prvr + -cjjr^ 

jabe jcde^A jabe jcde^^A 

+ P ^f + c A1 ^ v{p) + P ^y + c A a^] L ^ {p) (3 - 17) 

where we have introduced the subscript to avoid confusion with the set (|2.10p . with the propa- 
gator for p ab being unchanged. More specifically, in the Landau gauge the full set of propagators 
would be 



^acgbd jabe jcde jabe jcde^i 

(^ b (p)C d (-p))eS = ~ —f-V»» + CAp2 Vv» + p2 [( p 2 )2 + ( g A7 4] P ^(P) 
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xacsjbd 

(p»\pK d (-p))e$ = «(P)^(-P))eff = " -jr-Vr ■ (3-18) 
Taking the colour adjoint projection of gives 

(/^^(p)/^e s (-p))e ff = p2[(p2)2 ^ 7 4] ^(p) (3-19) 

from which one could, in principle, recover the original Gribov propagator when the leading 
order term of (|3.2p is applied to this. So one would have a hidden gluon with the perturbative 
propagator emerging as usual as 7 — > 0. In exploring this naive elimination idea further in 
general terms, and ignoring contributions from the path integral measure, it ought to be the 
case that in the infrared there is enhancement of various colour channels of the propagator 
as has been observed recently, [Ml EH]. If this were the case then (|3. 19j) might produce a 
hidden gluon propagator which freezes to a non-zero value as has been observed on the lattice 
by various authors, |25l|2Zl|2Sl|221l3niEIlE2ll33lEllE51- However, in speculating about the 
potential of a notional effective Lagrangian involving only the fields of (|2.6p which are enhanced 
in the infrared we must be clear in stating that such a theory has not been constructed. Our 
naive elimination by equations of motion, despite the fact that here it is related to the horizon 
constraint definition, is not the correct or accepted normal procedure to produce an effective 
quantum field theory. We offer it as a possible line of future investigation. In other words it 
might seem a natural way to proceed to study infrared properties by focusing on the actual 
fields which dominate at low energy. In toying with this idea it is certainly elementary to see 
that the Faddeev-Popov ghost, p ab and still remain enhanced at one loop. This is primarily 
because their associated vertices are effectively unchanged at this order. The difficulty comes in 
the sector where, although there is now an infinite number of interactions, which is not really 
a calculational obstruction, it is not clear whether renormalizability is retained. Therefore, one 
cannot even begin to consider if various colour channels of the ^ propagator enhance. However, 
in models of the low momentum the latter property is not crucial since, for example, the Liischer 
term is derived from a non-renormalizable construction, |46[I47|. Though an effective low energy 
theory with an infinite number of couplings of the colour valued fields to the quarks could 
be construed as a mimic of a flux tube model. Finally, irrespective of whether this naive use of 
an equation of motion is correct or not in trying to construct a theory involving only the fields 
which enhance in the infrared, it seems clear that one has to retain the horizon condition in 
some form and hence the gap equation for 7. The equivalences of (13.150 would appear to be a 
useful observation in this respect as we now have a vacuum expectation value definition of the 
Gribov horizon in terms of local function of only albeit an infinite perturbative series one. 
Though there may be a non-perturbative definition and thence an alternative way of pursuing 
an effective theory of enhanced fields. 



4 Static potential. 

In this section we concentrate on computing the static potential of heavy quarks in the three 
dimensional Gribov- Zwanziger Lagrangian. First, the calculational formalism was developed 
primarily for four dimensional QCD with massless quarks in the case of the canonical gluon 
propagator, [48l |49l [50]. It is based on the Wilson loop and a series of Feynman rules were 
constructed in coordinate space for the extreme case of a temporally long but spatially thin 
loop. With the one loop static potential emerging in [28l 091 [50], the two loop MS expression 
was produced later. This was derived first in the Feynman gauge, |51} [52], and then repeated 
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in an arbitrary linear covariant gauge in [531 154j. which verified the gauge independence of the 
potential. The latter computation is comprehensively detailed in |53| to which we refer the 
interested reader for background to technical aspects which we assume here. More recently, 
the three loop MS potential has been constructed in [551 E3 [571 EH E3 EO] which represents 
the current state of the art. By contrast only the one loop potential has been determined in 
three dimensional QCD in |53[ I6T] for canonical gluons. More specifically the momentum space 
potential, V(p), is, [53], f6T) . 



QCD 
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7 1 
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(4.1) 



Interestingly if one performs the inverse Fourier transform 

d 2 k 



V(r) 



(2nf 



V(k) 



(4.2) 



then the first term of (|4.ip reproduces the usual two dimensional Coulomb potential but the one 
loop correction rises linearly. However, as noted in [53} I6T] it is not clear what occurs for the 
two loop correction. Given the dimensionality of the coupling constant it is not inconceivable 
that higher loop corrections could lead to higher powers of the spatial separation r. It is worth 
noting that in four dimensions a linearly rising potential would correspond to a dipole term 
in the momentum space potential. By contrast, in three dimensions one requires a behaviour 
of 0(l/(p 2 ) 3//2 ) for a linearly rising potential, [53]. Given that we are interested in examining 
a theory, (12. 6j) . which is believed to be confining since it is consistent with the Kugo-Ojima 
confinement criterion, |12[ [T3] , our aim here is to compute the 7 dependent extension of (|4,ip 
and examine its behaviour in the infrared limit. If a term of the form 0(l/(p 2 ) 3 / 2 ) emerged 
then that would correspond to the confining potential being preserved in the presence of the 
Gribov horizon. Though for comparison we note that in [25] the analogous term, which would 
be a dipole, did not occur in four dimensions. One feature which emerged for the leading order 
term there was that the presence of the width in the gluon propagator meant that the coordinate 
space potential crossed the axis. Although this is present in the accepted form of the potential 
as computed say using lattice regularization, the potential actually crossed the r = axis at an 
infinity of places corresponding to a Friedel type of potential with a set of quasi-stable vacua. 
The key point was that the width was essential for this. A model where the gluon solely has a 
mass, if one ignores briefly the contradiction with the non-abelian gauge principle, would lead 
to a Yukawa potential in coordinate space which is never positive in r > 0. 

As we are extending the derivation of the static potential for the three dimensional version 
of (12. 6p we briefly recall several of the key parts of the formalism which were discussed in more 
detail in [25]. First, the potential is defined in terms of the Wilson loop which has a large time 
separation in comparison with the radial distance r, [481 1491 [50], 



V(r) 



lim — In ( 
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As is known this is equivalent to the definition involving the path integral, |48 [ I49 [ [50]. 

1 trZ[J] 



V(r) 



lim 

T^oo iT trZ[0] 



where 



Z[J] 



VA^VipVipVcVc exp 



d 3 x 



(4.3) 

(4.4) 
(4.5) 
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and the source term corresponds to placing heavy quarks according to 
J a Jx) = gv^T a (x - ir) - 8® (x - §r' 



(4.6) 



Here = tj^q is a unit vector which projects out the time component of the gluon it couples 
to. We also define r = |r — r'|. The presence of the sources introduces additional Feynman rules 
which are not dependent on the spacetime dimension and are given in, for example, |49} 153] . 
Though we follow the more modern approach and perform our static potential computations 
in momentum space rather than directly in coordinate space. The former is connected via the 
inverse Fourier transform, (|4.2p . The extension of the formalism to the Gribov-Zwanziger case 
is to replace the Lagrangian (|2.ip in the path integral in (|4.4p by the Lagrangian (|2.6p whence 
the measure is extended to include the localizing fields, 



Z[J] 



VAaV^V^VcVcViVpVujVuj exp 



d'x (L GZ + 



J a ^A a 



(4.7) 



As the localizing fields are completely internal they do not couple to the heavy quark sources. 
Hence at leading order the only field which is exchanged is the gauge potential which means 
that its propagator essentially determines the static potential at this order. The localizing Bose 
ghost plays no role until loop corrections are included. So, for instance, 



V( P ) 
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and performing the inverse Fourier transform gives 
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where Yq(z) is the Bessel function of the second type or Neumann function which is an entire 
function. Interestingly the four roots of the algebraic equation 



C A J 4 



(4.10) 



emerge as the arguments of the functions, 
behaviour noted in [531 161] 



Taking the 7 — > limit recovers the Coulomb 



lim V(r) 

7-s-O v ' 



CfI 
2tt 



ln(r) + 0(g 4 



(4.11) 



However, if one plots the functions of (|4.9p for non-zero 7 the coordinate space potential has a 
similar feature to the Friedel form of four dimensions. Although the Coulomb potential crosses 
the axis once when 7 is non-zero the static potential has an infinite number of crossing points 
which would again lead to quasi-stable vacua. 

We now turn to the one loop corrections of (|4.8p , A representative set of topologies for this 
calculation is given in Figure 1 where we have displayed the full set of corrections to the 2-point 
functions. In those first four graphs the blob represents all the one loop corrections. However, 
due to the mixed propagator the correction to the 2-point function occurs even though there 
is no direct source ^ coupling. The two box diagrams are important for the exponentiation 
implied in the definition of the static potential, (|4.4p . and the asscociated issues with have been 
discussed at length in |48 } 149 } [50 | 162 } 163] . Therefore, it remains merely to compute the diagrams 
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Figure 1: One loop topologies contributing to the static potential. 

explicitly. Essential to this is the automatic Feynman diagram package Qgraf, [37], where the 
graphs are generated electronically and then converted into Form input notation. The algorithm 
we use is to break up all the Feynman graphs into simple master integrals and then identify the 
explicit functions for three dimensional spacetime. A comprehensive analysis of such masters to 
three loops are given in [39]. However, for our situation there are only two main integrals but 
there is the complication of having to work with a non-standard propagator which induces the 
canonical propagator to have a width after application of simple partial fractions, such as 



where now the momentum can involve the external momentum. In four dimensions the explicit 
expressions for Feynman integrals involved masses, m, and momenta, p, appearing in the form 
m 2 and p 2 respectively. However, in our three dimensional case with the drop in one unit of the 
dimensionality of the integral measure, the dependence is purely in terms of m and \pp 2 . As 
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\p 2 + i^Cll 2 \ 2 [p2 - iVCU7 2 ] 



(4.12) 
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discussed earlier for fields with a canonical mass term this is not a significant issue but when 
there is a width present one has to find the square root of the corresponding squared mass 
of the propagator. We follow the procedure used previously but allowing for the presence of 
the external momentum. Again ultimately one should obtain real and not complex expressions 
which is a check on our reasoning. So, for example, we have used the following intermediate 
expressions 
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(4.13) 



by adapting the results of [39] in the same way as before, where J k = d k/(2n) . In (|4.13p we 
will use our simplification identities which allow us to write the functions of a complex variable 
in terms of a real and imaginary part. Although the three dimensional theory is finite, we still 
work in dimensional regularization with d = 3 — 2e as some of the master diagrams have poles 
in e, [39]. However, whilst the theory is ultraviolet finite we have been careful to check that 
there no (spurious) infrared infinities arise as a consequence of breaking the Feynman graphs 
up into scalar master integrals. For instance, such divergences could arise from the 1/p 2 part of 
the propagators in the transverse projector or its powers but again we have checked that such 
potential terms cancel among themselves. This finiteness, at least to one loop, ensures that 
there is no source gluon renormalization constant as there is in the four dimensional arbitrary 
gauge calculation. 



Given these considerations we are now in a position to record the one loop correction to (14.8 
for (ESI). We find 
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where we have defined the intermediate functions rji(p 2 ) by 
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Essentially these arise from taking the real and imaginary parts of expressions such as those 
given in f|4. 13j) . 



As a check on (|4.14p if we take the 7 — > limit we recover the one loop expression of |53t [6l~] 
given in (|4.ip . We stress though that whilst (|4.ip was computed in an arbitrary linear covariant 
gauge we have been restricted to the Landau gauge as we have incorporated the Gribov problem 
within the Lagrangian. However, ()4.14p represents the first non-trivial check on ()4.ip . Given 
that (|2.6p is supposed to represent a confining theory we can now examine (|4.14p to see if the 
dominant behaviour in the zero momentum limit could deliver a behaviour leading to a linearly 
rising potential. In three dimensions this would correspond to an 0(\fp 2 / (p 2 ) 2 ) type term and in 
(|4.14|) there is one term which appears with such a singularity. However, its numerator involves 
the function i]2(p 2 ) which vanishes at zero momentum and therefore, the appropriate singular 
behaviour does not seem to emerge. More concretely as p 2 — > we have 
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Therefore, similar to the four dimensional case, [25], the one loop correction freezes to a finite 
value, which is 



V(0) 



+ 0(g b ) 



(4.17) 



487T 7 3 

and there is no net divergence whose presence would at least be necessary for a rising potential. 
There is a degree of irony with this observation in that the 7 = potential at one loop has a 
Fourier transform which produces a linearly rising potential. Though in that case in is not clear 
what would transpire at two loops given the dimensionality of the coupling constant. However, 
as noted in [2H |25] the more appropriate route to proceed down would be to analyse the zero 
momentum behaviour of the propagators of the localizing fields. Whilst the fermionic ghosts 
are both enhanced at zero momentum in (12. 6p . it has recently been shown that the same is true 
for certain colour components of the Bose ghost, [Ml [25]. This property is independent of the 
spacetime dimension, [25] , but it has not been fully determined what the implications are for 
the static potential. 
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5 Formal propagator corrections. 



We turn now to the enhancement of the Bose ghost fields. Recently the structure of these fields 
was analysed non-perturbatively by Zwanziger in [23] where it was demonstrated that there 
was enhancement in certain colour channels. The result is based on the spontaneous breaking 
of the BRST symmetry but given the presence of the horizon condition, which equates to a 
constraint on the gluon and fields, this requires a more careful analysis than usual. One key 
outcome is that the associated Goldstone bosons of this spontaneous breaking generate massless 
excitations non-perturbatively but crucially in the context of the Gribov-Zwanziger Lagrangian, 
these fields are enhanced in the infrared. One consequence is that this enhancement is present 
order by order in perturbation theory and this was confirmed by explicit one loop calculations 
in four dimensions in the MS scheme, [25]. However, in [25] the main emphasis was on the 
transverse part of the propagators since the interest was in examining the implications for the 
static potential for heavy coloured objects. The longitudinal piece of the exchanged particle does 
not contribute in the particular configuration considered for the Wilson loop. In this section 
we provide the formal construction of the full propagator for the gluon and ^ fields prior to 
considering the behaviour in the infrared limit when the gap equation is realised which is given 
specifically in the next section. Although we concentrate primarily on three dimensions in this 
article we will also include the analysis for four dimensions in this section. This is because the 
general reasoning of |24] is dimension independent and we confirm this at one loop by considering 
both dimensions within our analysis here. 

We begin by recalling how enhancement occurs for the simple situation of the Faddeev-Popov 
ghost, c a . As was demonstrated in Gribov's seminal contribution, [5j, one computes the ghost 
2-point function at one loop when the horizon condition is implemented in the path integral. 
The consequent presence of the Gribov mass in the gluon propagator produces an expression 
which differs from what one would obtain using a canonical propagator. Expanding the finite 
function of the momentum in a Taylor series around zero momentum it transpires that the 
leading term is equivalent to the one loop correction to the Gribov mass gap equation. As the 
theory has no meaning as a gauge theory unless 7 satisfies the gap equation, this implies that 
the leading term of the 2-point function in this limit is not 0(p 2 ) but O ((p 2 ) 2 )- Hence, the low 
energy behaviour of Faddeev-Popov ghost propagator is not the usual perturbative form but the 
enhanced dipole form. Subsequently, this property of ghost enhancement has been translated 
into the Kugo-Ojima confinement condition. This was originally established in [12\ [T3] for 
the Landau gauge version of Yang-Mills involving only Faddeev-Popov ghosts. More recently 
the criterion has been extended to the Gribov-Zwanziger context where there are additional 
fermionic ghosts, which are clearly absent in the original Kugo-Ojima BRST analysis, as 
well as the Bose ghosts, p ab and [13] • In the context of (|2.6p it has now been established 
that iu® b also enhances. It has been demonstrated in the full analysis of [241 [36] and in explicit 
two loop MS computations in four dimensions, [23]. We note at this point that we have repeated 
the latter calculations for w^ 6 in the three dimensional version of (|2.6p and have verified that 
when 7 satisfies the two loop gap equation, (13.8p . the fermionic ghost enhances. 

In focusing on the Faddeev-Popov enhancement derivation the key ingredient is the zero 
momentum limit of the 2-point functions. As we have evaluated the 2-point function corrections 
exactly at one loop for all the fields we can now consider the situation for the Bose ghost fields. 
For ease we consider the p ab field first. In the zero momentum limit we have, 
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whence it is elementary to observe that the leading one loop correction is precisely that which 
appears in the gap equation, (|3.8j) . Hence, when 7 fulfils that condition the leading term of the 
2-point function is the O ((p 2 ) 2 ) part of the one loop term. Consequently when one inverts the 
coefficient of the tensor in this zero momentum limit then the p ab field enhances. One could, 
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Although the explicit expressions for the two loop corrections to the 2-point functions are not 
available, it is possible to determine their zero momentum behaviour using the vacuum bubble 
expansion. A similar procedure was followed in the four dimensional case and we note that at 
two loops the leading momentum part of the p ab 2-point function is precisely the two loop gap 
equation f|3.8|) . Therefore, the enhancement of p ab is present at next order in exact agreement 
with the general BRST arguments of [24J. 

Whilst this demonstrates that a part of the Bose ghost can enhance we need to complete 
the analysis by considering the imaginary component which is more involved due to it being 
entwined with A°^. In three dimensions the zero momentum limit of the 2-point function is 
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Clearly the piece analogous to (|5.1|) would equate to enhancement if one could perform the zero 
momentum inversion of (|5.3p in the absence of the additional colour channels. However, not 
only would this not be correct it would ignore the fact that the construction of the propagators 
actually requires to be included due to the mixed quadratic term of (|2.6p . In order to do this 
we consider the problem of the inversion in a general context. First, we recall the procedure for 
the transverse sector, [25], and define the matrix of colour amplitudes of the 2-point functions 

tab 



of the set of fields pf] by 





( 


X5 ac 


Uf acd 





jy{ab\cd} 




Uf cab 


Qabcd 







\ 








(-\abcd 



(5.4) 



where 



Qabcd 
Qabcd 



Q i5 ac 5 bd + ^jacejoae + ^jaoe jcae + g^ca + p^aogca + f^aap 

Q p 6 ac 6 bd + w p f ace f bde + K p f abe f cde + S p df cd + V p 5 ab 5 cd + T p 5 ad 5 bc 



pace fbde 



cabe rede 



•abed 



sab sred 



srad vbc 



(5.5) 
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and, 

d abcd = iTr^TfTXtf) (5.6) 

is totally symmetric. The subscript on the group generator T\ indicates that it is in the adjoint 
representation. We have omitted the common transverse projector from (|5.4p . Including p ab 
in the basis here may not appear to be necessary since (|5.4p is block diagonal, and we have 
treated it already, but it is relevant when we examine the longitudinal sector since not all the 
corresponding zero entries of (|5.4|) remain zero. In introducing general amplitudes for the colour 
decomposition we note that each represents the leading term of the 2-point function as well as 
the loop corrections. In each of the rank four colour decompositions we have included structures 
which do not arise in the explicit computations at one loop. Aside from ensuring we work with 
a complete basis such structures may occur at a higher loop order but they will also give us an 
insight into the effect such pieces have on the form of the inverse of the matrix which is the 
propagators. For this we formally define the inverse in a similar way with 



U {cd\pq} = Bf pcd V cd P q q £ g _ 7 j 




where 

v cdpq = V ^cp 5 dq + j^cpejdqe + ^jcde jpqe + ^cdpq + M ^cdgpq + tf^cqgdp 
v cdpq = Vp5 cp 5 dq + Jpf cpe f dqe + ^fcdejpqe + ^^dpq + J^^cdgpq + ^gcqgdp 



(5.8) 

and the transverse projector is again omitted. As n^"^ is the inverse colour matrix, it satisfies 

/ 5 cp 

A {ab\cd} u {cd\pq} = I q § cp 6 dq Q | (5.9) 

5 c P5 dc t 



\ 



where the right hand side is the unit matrix in the colour vector space of the basis of fields we 
use. For the inversion relating to the Lorentz structure, we recall the trivial identity 

Vpu = P^u{p) + V(P) ( 5 - 10 ) 
and note that the first term on the right side is where our current focus is. 

The method to find the formal inverse is to multiply out the matrices and solve the resulting 
relations between the amplitudes algebraically. In order to do this we note that the products 
involving d°^ cd can be simplified with the relations, [25], 

abed 



d abpq d edpq = ^gabged + ^ fgacgbd + gadgbc\ + fflg (jacejbde + f ade f bce^ + ^al 
f apejbqe d edpq = ^gabged + ^ (gaegbd + gadgbc^j 

+ b 3 (f ace f bde + fdefbcfj + hd abed ^ n 



where the coefficients are defined by, [25J, 

iabcd\ ^ 



0-1 



540C%N A (N A - 3)d$ cd d c * pq d £ pq + lte(2N A + 19) (df cd df c j 
- 150C\N A (3N A + ll)df cd df cd + 625C|iV"j 
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02 



54N A (N A - 3)[12(N A + 2)d e J 9h (fJ 9h - 25C\N A ] 
144(11 Aa - 8) (d a A bcd d a A bcd ) 2 - W80C 2 A N A (N A - 3)d a A bcd d c A dpq d a A bpq 



+ 625C A Nl 



abed mbed 



mmc\N A dT d a A 
1 



1<0&N A {N A - 3)[\2{N A + 2)d e [ 9h d A Sgh - 2hC\N A ] 



and 



as 



a 4 



2b 2 



[12(N A + 2)d a A bcd d a A bcd - 2hC\N A ] 
MC A N A (N A - 3) 

[21Q{N A + 2)d a A bcd d c £ pq d a A bpq 



125C A N A - 360C 2 A d A bcd d abcdl 



18[12(iV A + 2)d e [ gh d e l 9h - 2hC\N A ] 



(5.12) 



[5C\N A - 12d a A bcd d a A bcd ] 



bs = 



[6(Na 



l)d a A bcd d a A bcd 



5C A N A ] 



9C 2 A N A (N A - 3) 



9C A N A (N A - 3) 
Multiplying out the matrices results in the relations 

1 = AX + C A UB , = XB + (vz-Mz + CAKz + ^CAJijU 

= AU + (q 5 + + ^C A W ( - T^j B 

1 = + & 2 £?V^ + b 2 S^ + a 2 S ? £ ? + 
= 











b 4 = 



3 

(5.13) 



Qi + CaW 5 + -C%S i + NaV^ +Tt)Mt + C A ^ + + Z> € + m 



6 2 £ ? W C + 6 2 <%^ + a 2 S^ + Q ? W 5 + T 5 P f 

+ + S C Z> C + b 4 W^ + 64^^ + a 4 S^ 



V^P S + Q(Jt + -CUW^ + 2b 3 S^ + 2a 3 <%£ 5 + 26 3 W f £ e + + T 5 J" f 



2 ' ' 2 

b 3 S^ - a 3 S^C^ - T^J^ 



- Wt.Ni. + KtPi ~ 
for the and ^ sector. The full set for the p ab sector is 
1 = Q P V P + b 2 CpW p + b 2 S p J p + a 2 S p C p + T P M P 

= [q p + CaW p + ^C 2 a S p + N a V p + T^M p + (cAJp + ^ClCp + Vp + N^V, 

+ hW p C p + InSpJp + cnS p C p 
= b 2 CpW p + b 2 S p Jp + a 2 S p Cp + Q P W P + T p V p 

= Q P £ P + Wpjp + SpVp + b 4 W p Cp + b 4 S p Jp + a 4 S p £ p + SpM p + T p £ p 
= WpVp + Q pi 7 p + X -C A W p J p + 26 3 S pi 7 p 
+ 2a 3 5 p £ p + 2b 3 W p C p + W p Afp + T P J P 



(5.14) 



Q P /C P + ^CaW pi 7 p + X -C A WpKp + ^Ca^pJp + C A lZ p lCp 
- W P AA P + 7e p P p - ft p A/" p - 6 3 5 P J" P - a 3 S p £ p - T P J P - 



- b 3 W p Cp 

TplCp . (5.15) 
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Whilst this set appears formally similar to the corresponding equations of the A® and subset, 
the final equation has one fewer term. 

For an arbitrary colour group the full solution to both sets of equations are cumbersome. For 
instance, by way of illustration for the slightly simpler case of the sector we have recorded the 
full expressions for SU(N C ) in Appendix C. For the top sector the gluon and mixed propagators 
are simple for an arbitrary colour group giving 



A 



B = - 



[gg + C A Kt + |gAVVg] 
[(Q ? + CaR* + \C A W^)X - C A W] 
U 



[(Q ? + CaK^ + \C A WAX - C A U 2 \ ' 
However, for the propagator we only present the expressions for 577(3) which are 
1 



(5.16) 



2Q 6 



3(3,% - 2Wz)(St + 2W 5 )(<% + Wf) + 8(7<% + 3W ? )Q| 
+ 16Q| + 2(275| + 205 c >V € -8>V|)Qf] 



ICp = 



4Wg 

~ [2Q ? + 3c% + 3W 5 ][2Q ? + 3S f - 2W 5 ] 

^- [(4(3<% - Wz)Qt + 3(3<% - 2W ( )(S, + W^))(2U 2 - W,X - 27%X) 



Mi 



- 8(KtX -U) 2 Q\ 

x 2Q,X + 6TZpX -GU 2 + 3W^x] _± [2Q f + 3<% + 3W ? ] _1 [2Q 5 + 3<% - 2Vty _1 
- 4 [2Q& + (3<% + 2V^)(«% - V^)] 
x [2Q f + 3<% + SW^ 1 [2Q ? + 35 ? - 2W 5 ]~ 1 [2Q 5 + S ( + 2W c p 1 

6 [215| + 5|W{ - 12<% Wf - 4W| + 2(7<% + 4W c )Q e <%] [2Q f + 15<% + 6W ? ] _1 
x [2Q C + 3<% + 3W5]" 1 [2Q 5 + 3<% - 2W 5 p 1 [2Q C + + 2W € ] -1 
[(3<% - 2V^)(5f + 2W ? ) (<% + W ? ) + 2(<% + 4W C ) Q^] 



2Q ? 



[2Q C + 3c% + 3^]^ [2Q 5 + 35 c - 2W^]" 1 [2Q C + + 2W ? p 1 



(5.17) 



Next we formally repeat the procedure for the longitudinal part of the matrix of 2-point 
functions and denote the corresponding quantities with a superscript L . For instance, the matrix 
of 2-point functions is now 



where 



QL abed 
L abed 



Q 



p^L {ab\cd} 



f X L 5 ac 
U L f cab 
{ V L f 



)L facd 



■j^L jacd ja 
qL abed q 
L fcab q r\L abed 



(5.18) 



QLgacgbd + y^L jace jbde + jabe jede + g^abed + pL^ab^cd + j-Lgadgbc 
QL^ac^bd _|_ yyi jace jbde _|_ j^L jabe jede _|_ gL^abcd _|_ ^pL^ab^cd _|_ q-L ^ad ^bc 

(5.19) 
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In writing the formal longitudinal part we are making no assumptions at the outset concerning 
the form of the 2-point functions. For instance, from (|5.ip and (|5.3p it is clear that 



Qi = Q\ + 0(a 2 ) = Q p + 0(a 2 ) = Q L p + 0(a 2 ) (5.20) 

to one loop but we do not impose that condition initially. There is a non-zero entry in the A^-p 1 ^ 
slot due to a non-zero one loop contribution to the longitudinal part of this 2-point function. 
Thus the inverse has to be more general than that for the transverse sector and we take 

/ A L 5 cp B L f cpq C L f cpq \ 

Y[L{cd\pq} _ fiL pad j^Lcdpq gLcdpq j (5 21) 

y C L f pcd £^ cd PQ gLcdpq J 

where 

v Lcdpq = V L 6 cp 5 dq + jLjcpejdqe + jcde jpqe + j^L^dpq + J^Lgcdgpq + J^f^cq^p 

£ Lcdpq = £ L 6 cp g dq + jrLjcpejdqe + gLjcdejpqe + yLjfdpq + yLgcdfiM + Z^5 Cq 5 dp 
£ Lcdpq = £ L 5 cp 5 dq + jrLjcpejdqe + gL jcde jpqe + ^Ljcdpq + yLgcdpq + Z^cqgdp 

(5.22) 

The inverse containing the longitudinal sector of the propagators satisfies a similar equation to 
that for the transverse sector, 

/ 5 cp \ 

A L{ab\cd} u L{cd\pq} _ q ^cpfidq q _ (523) 

\ 5 cp 5 dq ) 

However, in order to reduce the size of the algebraic equations we will have to solve eventually, 
for the longitudinal sector we set 

W p L = K L p = 5 p L = V L p = T p L = (5.24) 

at the outset since it is evident from the explicit computations at one loop that these relations 
are valid. If at higher loop order it turns out that any of these is non-zero then one would have 
a different set of equations to solve for the propagators. We find 



1 = A L X L + C A U L B L + C A V L C L 

= x l b l + (vf -M[ + c a k\ + x -c A j^u L + Uf - z£ + c A g\ + \c A ^) V 
= x L c L + ^ - z\ + c A g\ + \ca?£\ u l + [s L p - z L p + c A g p + \c A ^ v L 

= A L U L + (q% + C A 1Z\ + X -C A W\ - rfj B L 

1 = Q\v\ + b 2 C\w\ + b 2^Jt + a 2 S^Cl + TfM[ 

= (q% + c a w^ + ^c 2 a s% + n a v% +tAm% + (c A Jz+\c 2 A c% + V^+M^jV^ 

= b 2 £fw£ + b 2 S^J^ + a 2 S^C\ + Q^M[ + T^V\ 

= Q\L\ + W^JjF + S\V\ + b^N\L\ + hSj;^ + a A S^Cl + S^Aff + T£l\ 

1 
6 
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+ 2a 3 S££% + 2& 3 W|x£ 



6 

,L r L 













2 — 2 



Q\ + CUWf + -C\S\ 

LnjL i l qL -nL i njL r>L 



+ c a k\q\ - b 3 w^nf - w\z\ + - n\z\ 



V L A L + Q^C 



LnL 







Qp£ L 



Q L P K 








V L B L + Q^Ql 
Q L P fp , o 



V L C L + Q L p G p 



■£p '-p 

iLr>L 








Q£4 
2^ 








Q L P z L P 

(5.25) 



These are clearly more involved than the transverse sector and lead to more complicated forms 
for the explicit amplitudes. In addition to our earlier nullifications, with T^ L = at the outset 
we find 



A 1 



Qp[(Qc + c A n L + \c A wk)x L - c A {u L y\ - c a [q l + c A iz\ + \cM\(y L y 



B = — 



U Q 



Q L P [{Q L t + c A n\ + ±c A nf )* L - c A (u L y\ - c a [q l + c A n\ + \c A w^]{v L f 

[Q L + c A n\ + \c A y^l\v L 



Qp[{Q\ + c A n\ + \c A w%)x L - c A (u L ) 2 } _ c a [q\ + c A n\ + \c A w%\{v L y 

(5.26) 

for an arbitrary colour group. For the remaining amplitudes restricting to SU(3) produces 



V? = 



1 



2Q\ 



3(35^ - 2W^){S^ + 2W ? L )(^ L + W 5 L ) + 8(75 ? L + 3Vtf )(Q[) 2 
+ 16(Qf ) 3 + 2(27(5 5 L ) 3 + 205 5 L W C L - 8(W f L ) 2 )Qf 



2Q^ + 35f + 3Wf [2Q^ + 35 f L - 2Wf [2Q^ + 5f + 2Wf 
4Wf 



-1 r 



[2Q^ + 35 f L + 3W f L ] [2Q^ + 35 f L - 2W f L ] 



^ = - 



1 

Of 



8(Qf) 2 Q^^ L - 8(Q^) Z Q^(^) Z - 24(Q^)^^(V^) Z + 2AQ%Q$n%S£X 



■^L\2 r\L n ,L\2 
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Ml 



+ mlQ L p ^{u L f - AQ^Q L p X L {W^) 2 - 72Qfnfsl'(y L ) 2 

+ 24Q^W^(V L ) 2 - 36Qf S^(V L ) 2 + \2Q\ (V L f(W%) 2 

+ \%q%k\(s%) 2 x l + §Q L p n\s\wlx L - l2Q L p nl{W^) 2 X L 

- \m L p {S^) 2 {U L ) 2 + 9Qj>Vf (Sf ) 2 * L - QQ L P W^{U L ) 2 S^ 

+ m L p s^{w^fx L + i2Q L p {w^) 2 {u L f - m L p {^fx L 

- 54^(5|') 2 (v L ) 2 - isn^{v L ) 2 w^ + mn^{v L ) 2 {w^) 2 

- 27W^{V L ) 2 {S^) 2 - 9S^(V L ) 2 {W^) 2 + 18(V L ) 2 (Wf ) 3 
2Q\Q L p X L - QQ L p {V L f + 6Q%K%X L - 6Q L (U L ) 2 + 3Q^W^X L 



- 4 



18^(V L ) 2 - 9(V^ 2 



4 Wf 



-1 r 



2Qf + 35f + 3VWf 



-1 r 



2Qf + 35 e L - 2W $ L 



2Q^«Sf + (&Sf + 2Wf )(«Sf - VWf 



2Q 



L 
?L\3 



3S^ 



3Wi 



-i 



-l 



2 Q\ + 35^ - 2 VWf 2 + 5f + 2 VWf 



21(S£f + ($f ) 2 Vt£ - 12S^(W^) 2 - 4(Vt^) 3 + 2(75^ + 4Vt£)G£S£ 



,£\2 



L\3 



2Qf + 155 5 L + 6W ? L 2Qf + 35f + 3W| 



-l 



2 + 35f - 2 VWf 2 + 5f + 2 W e 



-i 



2Qf 



(35^ - 2W^)(«S£ + 2V^)(«S^ + Vt^) + 2(S£ + 4W^)Qf 5 f L 



x 



2Q[ + 35^ + 3Wf 



2Qf 



35 f L 



2Wi 



-i 



2g 5 + c% + 2W 5 p 1 



5 ~ 
2W L V L 



ft 



2[Q^ - 3(V L ) 2 ]Q^ - 9[2ft£ + ] 



- 3[2(W Z ) 2 -W^X 1 * -ZRjfX^Q 



Q L 



ft 



y 



p '-p ~p 

>£ l •\A)L-\ i nni /■\)L\2/ / ^)L\— 1 



3[2ft£ + Wf] + 2Q^j (VY(Qp ) 
2[Q^ L - 3(V L ) 2 ]Q£ - 9[2ft£ + W ? L ] - 3[2(ZY L ) 2 - W%X L - 2U^X L ]Q^ 



(5.27) 



As a check on these solutions we have verified that the actual propagators, ()2.10j) . are correctly 
reproduced when the a independent values of the 2-point function are inserted. 



6 £® and enhancement. 

Equipped with these solutions for both sectors we can now examine the specific problem of 
enhancement. For the transverse sector it is a straightforward exercise to substitute the explicit 
zero momentum behaviour from the 2-point functions (|5.ip and (|5.3p . As was noted in [25] this 
produces an enhanced £ p b propagator in the transverse sector as expected given our parallel 
reasoning for the Faddeev-Popov ghost propagator. Moreover, the enhancement in the three 
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dimensional case is similar to that of the four dimensional analysis, [25]. For the longitudinal 
sector the derivation of the zero momentum behaviour for requires care. Whilst the explicit 
expressions for the colour channel amplitudes in the longitudinal sector are general they mask 
the fact that ultimately we are in the Landau gauge. As is well known in order to construct 
the Landau gauge propagators in the non-Gribov scenario the gauge fixing term includes a 
parameter, a. When this vanishes one is in the Landau gauge. However, such a term is required 
in order to prevent a non-singular matrix inversion such as that needed for deriving (12. 8ft . As a 
occurs in the longitudinal term in that case, we cannot omit it in the inversion for the full set 
of longitudinal 2-point functions. Specifically, the leading term of X L is the only place where a 
appears. However, in order to extract the correct zero momentum behaviour of the propagators 
one must be careful in taking the limit to the Landau gauge. It transpires that this must be 
taken first in all our expressions for the propagator amplitudes and then the zero momentum 
limit taken. The order of the limits is not commutative. Though to examine the problem of 
enhancement we must set the gap equation initially. Following this procedure we find the zero 
momentum behaviour of the propagators is 



(CWC(-p)) 



15^/2tt7 5 
30 v 7 ^ 5 



ac rbd 



30V27T 7 5 



f abe f Cde P,v(p) 



J A 



C 



V 4 ,' 2^2„2 



(6.1) 



Clearly the colour structures of the transverse and longitudinal parts of the propagator 
are different. However, if one contracts either field of either propagator with a structure func- 
tion then the enhancement disappears. This loss of enhancement for this colour projection is 
completely in accord with Zwanziger's all orders observations from the BRST symmetry consid- 
erations in [23]. Though it should be noted that there are 0(l/p 2 ) pieces which remain. More 
specifically, retaining the next term of the expansion as p 2 — > 0, we have 



(A^p)Al(-p)) 
(A a ,(p)Cl c (-p)) 

<£?(p)C(-p)> 



V2p 2 » 2 



9 



5 ab P,u{p) 



3847rCyV 



CaI 2 




1 + 



127T7 



f abc p^(p) 



15^2^ 

CT(P 2 ) 2 9 2 
~ 2\cTp 2 9 2 



6C A ai - \2C\b 3 + C\bi - 432C> 3 6 4 

- 72C\a 3 - 108a 2 6 3 + 9C A a 2 b 4 

- 9C A a 4 b 2 + 432C A a 4 b 3 + 108a 3 6 2 



x [12a 2 6 3 - CAa 2 b 4 - 12a 3 6 2 + CAa 4 b 2 

1 5V2TT-f 3 



-1 



8 ac 5 bd PAp) 



CT(P 2 ) 2 9 2 



5\/2tt7 
336C5' ' 4 P 2 g 2 



42Cja 4 - MC\h + 7C A b 4 - 3024C A a 3 b 4 
- 504C A a 3 - 1728a 2 6 3 + UAC A a 2 b A 
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144C A a 4 & 2 + 3024C A a 4 6 3 + 1728a 3 6 2 



x [12a 2 6 3 - Cao-iW - 12a 3 6 2 + C A a 4 b 2 ] 1 



8 ac 8 bd L lll/ {p) 



+ 



+ 



20^/27rCy 4 7[72a 3 - 6C A a 4 + 12C A fr 3 - C 2 A b 4 \ 
7[12o 2 6 3 — C A a2& 4 - 12a 3 6 2 + C A a 4 b 2 ]p 2 g 2 



r e f bde PAp) 



5^/27rC2 4 7[72a 3 - 6C A a 4 + l2C A b 3 - C\b 4 \ 
2[12a 2 6 3 - C A a 2 b 4 - 12a 3 b 2 + C A a 4 b 2 ]p 2 g 2 
30\/2tt 7 3 

CT(P 2 ) 2 9 2 

- 1 ° 7 ^ 7r7 [6C>4 - 12C A 6 3 + C A 5 4 - 3C A a 2 b 4 
7C A ' p 2 g 2 

+ 36a 2 6 3 — 72C A a 3 — 36a 3 6 2 + 3C A a 4 b 2 



+ 



x 12a 2 6 3 - C A a 2 b 4 - 12a 3 6 2 + C A a 4 b 2 

5V27rCy 4 7[72a 3 - 6CUa4 + 12C A fr 3 - C%& 4 ] 
2[12a 2 6 3 — CAa 2 &4 — 12a 3 6 2 + C A a 4 6 2 ]p 2 g 2 



f abe f cde p,Ap) 



f abe f cde L^(p) 



+ 



+ 



120^Cy 4 7 [C A b 4 - 12C A fr 3 ] 
7[12a 2 6 3 - C A a 2 6 4 - 12a 3 6 2 + C A a 4 b 2 ]p 2 g 2 

\^cTl\C A b,-\2C A b z \ bcd 
[12a 2 6 3 - C A a 2 b 4 - 12a 3 fe 2 + C A a 4 & 2 F<7 2 A ^ ' 



j_ 120\/27r7[12ai6 3 - C A a\b 4 - 12a 3 bi + C A a 4 fri] ^ab^cdp ^ 



+ 



7C 3 J 4 [12a 2 b 3 - C A a 2 b 4 - 12a 3 6 2 + C^fcJpV 
15\/27T7[12ai6 3 - C A ai& 4 - 12a 3 &i + C A a 4 b\] 6 d 



2<7 4 /4 [12a 2 6 3 - C A a 2 b 4 - 12a 3 6 2 + C A a 4 b 2 ]p 2 g 
15\/27r7 3 



+ 



5 ab 5 cd L, u (p) 



CT(P 2 ) 2 9 2 
5\/27r7 



+ 



21Cl /4 p2 5 2 



6C> 4 - 12C A 6 3 + C A 6 4 - 432C A a 3 6 4 

- 72C A a 3 + 108a 2 6 3 - 9C A a 2 6 4 
+ 9C A a 4 6 2 + 432C A a 4 6 3 - 108a 3 6 2 



+ 



x [12a 2 6 3 - C A a 2 b 4 - 12a 3 6 2 + C A a 4 b 2 ] 1 
15V27T7 3 



5 ad 5 bc P^{p) 



5\/27T7 



+ 



336C^V</ 2 



42C> 4 - 84C A 6 3 + 7C A 6 4 - 3024C> 3 6 4 

- 504C A a 3 + 1728a 2 5 3 - 144C A a 2 6 4 
+ 144C A a 4 6 2 + 3024C A a 4 6 3 - 1728a 3 6 2 



x [12a 2 6 3 - C A a 2 b 4 - 12a 3 6 2 + C A a 4 b 2 ] 1 



5 ad 5 bc L^(p) 



27 



(W(-p» 

(pt(v)p?{-v)) 



cd i 



30\/2tt7 3 

CT(P 2 ) 2 9 2 
30\/2Vy 3 



+ 



+ 



30\/2tt7 
30\/27r7 



(6.2) 



In order to compare with the colour adjoint projection of [23] the leading behaviour of the 
bosonic ghost in the zero momentum limit for an arbitrary colour group is 



r q f hrs {e«{p)c s {-p)) 



- 5 



•ab 



V u b3V2C 5 /g 2 _ [21C74 - 64rWV/]yW 



- 5 



ab 



(P 2 ) 2 9 2 

100V27T 



384 7 3 

' 3/4 - 3 30v/2CY 4 7r7 



5127 4 



+ 



7p 2 g 2 



+ 



AdC^g 2 ] 
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Not only does the enhancement disappear but the simple pole is also absent leaving a finite non- 
zero value at zero momentum. Although the tree term is absent at zero momentum the loop 
corrections lead to a finite value in this limit. Whilst it is tempting to assert that the freezing of 
the transverse part of this correlation of a spin-1 field carrying one adjoint colour label is what is 
observed on the lattice and regarded as the frozen gluon propagator of the decoupling solution, 
the absence of a transverse propagator would exclude this as an alternative explanation. As far 
as we are aware numerical work observes a tranverse infrared frozen gluon propagator with no 
longitudinal part which is enhanced or otherwise. Moreover, the enhancement of the Faddeev- 
Popov ghost still remains contrary to what the lattice observes, [26J [27l [28], [29] [30] EH [32] . Just 
for completeness if we perform the same colour contraction on the original propagator, (|2.10D . 
we have 
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So that the effect of implementing the gap equation in deriving the zero momentum behaviour of 
the propagator, appears to reduce the momentum structure of both the transverse and longitudi- 
nal components by one power of momentum for this specific limit. Though writing the original 
propagator with the adjoint projection in this way demonstrates the existence of a massless 
pole only in the longitudinal sector. So given that massless poles seem to lead to enhancement 
in other situations, (|6.3|) seems to be consistent with this observation. As an aside we draw 
attention to the contrasting structures of (|6.4p and our earlier naive effective propagator, (|3,19p . 

Whilst we have concentrated on the zero momentum behaviour for the propagators when 
the gap equation has been implemented, it is worth noting some general features of the full one 
loop corrections to the gluon propagator itself. In (|6.2p we gave the leading order behaviour 
of the gluon. Unlike the localizing fields there is no divergence in the zero momentum limit. 
Moreover, the propagator vanishes at one loop similar to the original propagator derived from 
the Lagrangian. This is in keeping with [24] which showed that the gluon form factor vanishes 
and hence is the key to showing positivity violation for ()2.6|) . That our calculations reproduced 
this at one loop is in fact a consistency check on |24j . However, given the form of the equations 
for the longitudinal sector of the previous sector, it also turns out that there is no longitudinal 
component for the gluon propagator which therefore remains transverse at one loop similar to the 
original propagator. In fact this is due to the longitudinal correction being proportional to the 



28 



gauge parameter which vanishes for the Landau gauge. Hence, these remarks are independent 
of the dimension and the same property is present for the four dimensional gluon propagator 
in (|2.6p . Equally the mixed A^-^ propagator remains transverse at one loop in the Feynman 
gauge in either spacetime dimension. 

As [25] only considered the enhancement of the transverse propagator, because the focus 
was on the implications for the static potential, we also record the situation with the enhanced 
propagators in four dimensions. At leading order in the zero momentum limit we have 
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Thus the propagators have the same colour structure as the three dimensional case and so the 
colour adjoint projected fields are clearly not enhanced in this situation either. Including the 
subsequent term in the series we have, for a general colour group, 
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For the colour group 577(3) the Bose ghost propagator at 0(l/p 2 ) becomes 
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as p — > 0. Repeating the analogous calculation to (|6.3p . we find a similar structure in four 
dimensions since 
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for the leading order behaviour for each structure. Again the transverse enhancement disappears 
and overall the transverse projection also freezes to a non-zero value. There is longitudinal 
enhancement in keeping with |24| and our observation on the massless poles in the original 
propagator. 



7 Discussion. 



The main motivation of the article was to provide the loop analysis of the three dimensional 
Gribov-Zwanziger Lagrangian, (j2.6j) . to the same order in perturbation theory which is currently 
available in four dimensions. In having achieved this we note that many of the key properties are 
preserved. For example, the enhancement of the fermionic ghosts and certain colour channels 
of the imaginary part of the Bose localizing field is evident. Moreover, we have provided the 
complete analysis of the construction of the latter to one loop for both the transverse and 
longitudinal parts. This additional enhancement of a Bose field is in keeping with the Kugo- 
Ojima ethos underlying (j2.6l) . [121 [131 EH G3] , which has been argued to be a necessary criterion 
for confinement. However, the actual mechanism of how this is realised in practical terms in (|2.6p 
is still not resolved. In the pioneering ideas of Mandelstam and others, [651 166] 16?] 168]. for four 
dimensional Yang-Mills it was believed that the rising potential was due to the single exchange 
of a colour valued field between heavy quarks whose propagator was of a dipole form in the 
infrared. Then it was the matter of a simple Fourier transform to coordinate space to produce 
the rising potential. With the emergence of the enhancement of certain colour channels of the 
propagator observed in [231 E5] the exchange of this colour quanta was explored to see if 
the dipole behaviour emerged. However, it transpired that the vertex colour structure nullified 
the colour structure of the associated dipole term of the infrared part of the propagator. In 
addition the momentum dependence of the vertex would also have led to a diminishing of the 
power of the momentum in the exchange. In repeating the analogous analysis here we have 
confirmed at one loop an underlying feature of Zwanziger's all orders BRST reasoning. That is 
the enhancement of the Bose ghost is independent of the spacetime dimension. In other words 
one retains the dipole behaviour in the infrared. If one were to have the exchange of a dipole 
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type term to produce a linearly rising potential in three dimensions then it is clear that the 
single Bose ghost exchange would not be the simple explanation. This is simply because one has 
to have a l/(j> 2 ) 3//2 behaviour in the infrared for the zero momentum limit in order to have a 
linear dependence on the radial distance upon performing the Fourier transform. Of course, this 
is on the understanding that the underlying mechanism is effectively dimension independent. 
To manufacture the necessary extra powers of momentum to alter the enhanced form via say 
vertex correction momentum dependence would appear to be difficult because to be balanced 
one would have to have two factors of (p 2 ) 1//4 - 

Whilst anomalous dimensions of vertices and fields can in principle acquire large corrections 
in the infrared, to explore this further would require a summation of a significant set of Feynman 
diagrams, for instance. Moreover, this would seem an unlikely avenue since the three dimen- 
sional theory is superrenormalizable being less ultraviolet divergent than the four dimensional 
counterpart. So an anomalous dimension could be trivial in three dimensions. Instead the point 
of view might be that the actual dynamics of how the rising potential emerges in (|2.6p rests in 
the enhancement of the Bose ghosts residing within Feynman diagrams themselves. Several ways 
to perhaps achieve this could be worth considering. One might be the summation of ladder type 
diagrams which could be similar to a colour flux tube. An advantage of this is that the dimen- 
sionality of the exchange required in both dimensions could be naturally accommodated by the 
Feynman integral measure. We have already touched on another possibility when we rewrote the 
definition of the horizon condition purely in terms of the Bose ghost which enhances. One could 
conceive of some sort of effective infrared theory involving only the fields which enhance. The 
canonical gauge potential which determines the ultraviolet dynamics would appear as a bound 
state of the Bose ghost, . Such a scenario of the gluon or gauge potential being regarded as 
a bound state has been considered before but in other contexts. For instance, in [69] the gluon 
was interpreted as the bound state of quarks. Whilst the analogous interpretation here would 
be a bound state of fields, the work of [69] demonstrated that the (perturbative) structure 
of d-dimensional QCD was equivalent to the non-abelian Thirring model for calculations. Al- 
though this was primarily true only at a non-trivial fixed point of the renormalization group 
flow in d-dimensions, one could in effect use the simpler non-abelian Thirring model. In other 
words it was an effective theory but which is non-renormalizable in dimensions greater than 
two. It is not inconceivable therefore, that underlying the Gribov-Zwanziger theory there is a 
similar but clearly more complicated non-renormalizable effective theory which could have the 
structure of a nonlinear a model with an infinite set of interactions. In essence the quark gluon 
vertex could be replaced by quarks interacting with n fields. However, it is difficult to see 
how such speculative ideas could be realised practically in the short-term. In passing we note 
the curiosity that the non-abelian Thirring model has a formulation which involves a dimension 
two operator built from a spin-1 auxiliary field carrying an adjoint colour index. 
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A Transverse parts. 



In this appendix we record the explicit one loop contributions to the transverse parts of the 
2-point functions. We have 
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B Longitudinal parts. 

In this appendix we record the explicit one loop contributions to the longitudinal parts of the 
2-nnint functions. We have 
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Analogously to the previous appendix, we record the respective zero momentum limits are 
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C propagator. 

In this appendix we record the formal solution of (|5,15p for the case of SU(N C ) as those for an 
arbitrary colour group are too involved. Unlike the solution given earlier for the A p and £ p b 
sector we have not set T p to zero at the outset. We have, for the transverse sector only, 
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+ 72iV 2 7e p 5 p (T p ) 2 + 72iV 2 7e p r p (W p ) 2 + 36iV 2 (5 p ) 3 W p + 84iV 2 (S p ) 2 T p W p 

+ 24iV 2 5 p (T p ) 2 W p + 36iV 2 T p (W p ) 3 - 36iV c (Q p ) 2 (5 p ) 2 - 72iV c (Q p ) 2 (W p ) 2 

+ 57W C Q P K P S P W P + 432N c Q p K p T p W p + 72iV c Q p (5 p ) 2 T p + 288iV c Q p cS p (W p ) 2 

+ 360N C Q P T P (W P ) 2 + 216N C TZ P (S P ) 2 W P + 576N C TZ P S P T P W P + 432N C TZ P (T P ) 2 W P 

- 216N c K p (W p f - 36N C (S P ) 2 (T P ) 2 + 108iV c (5 p ) 2 (W p ) 2 + 288iV c cS p T P (W p ) 2 
+ 144iV c (T p ) 2 (W p ) 2 - 1087V C (W P ) 4 + 432(Q p ) 3 ft p + 864(Q p ) 2 7e p 5 p 



+ 1296(Q p ) 2 ft p T p + 432(Q p ) 2 S p W p + 864(Q p ) 2 T p W p + 432Q p ft p (S p ) 2 
+ 1728Q p ft p S p T P + 1296Q P ^ P (T P ) 2 - 432Q P ^ P (W P ) 2 + 216Q P (S P ) 2 W P 
+ 864Q p S p T p W p + 864Q p (T p ) 2 W p - 216Q P (W P ) 3 + 4327e p (S p ) 2 T p 
+ 864ft p S p (T p ) 2 + 432ft p (T p ) 3 - 432ft p T p (W p ) 2 
+ 216(5 P ) 2 T P W P + 432S p (T p ) 2 W p - 216T P (W P ) 3 

6Q P + N 2 S p + 3iV c W p + 6T P ] ~* [2N C K P + N C W P + 2Q p - 27^] _1 
x [QQ P + (N c + 6)5 P - 6W P + T P Y l [(N c - 6)5 P - 6Q P - 6T P - 6W P ] _1 [Q P - Tp]' 1 
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N 2 (S P ) 2 - 3N C S P W P + 18Q P 5 P + 18(5 P ) 2 + 18S P T P - 18(W P ) 2 



6Q P + N 2 S P + 3N C W P + 6% 



n -1 



x [(JV C - 6)5 p 

-6iV 4 P p (5 p ) 2 - 5iV 4 (5 p ) 3 + 18iV 3 7V> p W p + 9N%(S, V2 



168N 2 V p (S^ 2 



[6Q P + {N c + 6)5 p - 6W P + T,}- 1 

?;Q p -6r p -6W p ]- 1 

~ ' P S P W P + 9iV 3 (5 p ) 2 W p - 144iV 2 P p Q p cS p 

pj - UAN 2 V P S P T P + 108iV 2 P p (W p ) 2 - 126iV 2 Q p (5 p ) 2 
- 1447V 2 (S p f - 126iV 2 (S p ) 2 T p + 108iV 2 5 p (W p ) 2 - 108iV c P p Q p W p 
+ 72iV c P p 5 p W p - 108iV c P p T P W p - 216N C Q P S P W P - 108iV c (5 p ) 2 W p 
- 21QN C S P T P W P + 108iV c (W p ) 3 + 216P P (Q P ) 2 + 432P p Q p S p + 432P p Q p T p 
+ 216P P (5 P ) 2 + 432P p 5 p T p + 216P P (T P ) 2 - 216P P (W P ) 21 

6Q P + N 2 S P + 3iV c W p + 6T P ] ^ 

6(AT 2 - 6)V p + 5iV 2 5 p + 6iV c W p + 6Q P + 6T P 1 



x [6Q P + (N c + 6)cS p - 6W P + T,}- 1 [(N c - 6)S P - 6Q P - 6T P - 6W, 



i-i 
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K = \ {-N?(S P ) 3 + 9N*(S P ) 2 W P + 18iV 2 Q p (S p ) 2 + 72N 2 Q p S p T p + 36iV 2 (S ( 



x 



+ 1UN 2 (S P ) 2 T P + 72N 2 S P (T P ) 2 + 216N c Q p S p W p + 216N c Q p T p W p 

+ 108N C (S P ) 2 W P + 360N C S P T P W P + 216N C (T P ) 2 W P - 108N c (W p f 

+ 432(Q p ) 2 T p + 864Q p S p T p + 864Q p (T p ) 2 + 432(5 P ) 2 T P 

+ 8645 p (7;) 2 + 432(7;) 3 - 432r p (W p ) 2 " 

'6Q P + N 2 S P + 3iV c W p + 6T P ] ~* [6Q P + (N c + 6)S P - 6W P + T^ 1 
x [(iV c - 6)5 P - 6Q P - 6T P - 6W P ] _1 [Q p - r p ]^ x . (C.l) 

One of the reasons for recording these expressions is that they illustrate how cumbersome the 
full propagator structure for the two colour index field is even in the case where there is no 
matrix of 2-point functions. However, a more important feature is to illustrate the role of 
the implementation of the gap equation to extract the enhanced propagators. In these final 
expressions the enhancement will derive from the factors [Q p — 7^] _1 in the T> p , fC p and J\f p 
channels. If it transpired from explicit calculations at any loop order that T p was non-zero or 
not proportional to the gap equation itself then there would be no enhancement of p a p . The same 
situation would occur for but in respect of 7g in that case independent of the complication 
of having to analyse the 2x2 matrix for the gluon and £ ab sector. Similar comments apply to 
the longitudinal situation as the equations for the longitudinal piece are formally similar. 
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